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Radio timing observations of the high-magnetic-field pulsar PSR J1718— 3718 
have shown that it suffered a large glitch with Aug/u = (33.25 ± 0.01) x 10~^ 
between 2007 September (MJD 54336) and 2009 January (MJD 54855). This 
. is the largest pulsar glitch ever observed. As is common, there was a small 

. increase in braking torque at the time of the glitch but, unlike all other pulsars, 

. the braking torque has continued to increase over the two years since the glitch. 

Polarization observations show that the mean pulse profile has about 30% linear 
2 ■ polarization with a smooth change of position angle through the pulse and give a 

. rotation measure of —160 ±22 rad m~^. There was no detectable change in pulse 

profile at the time of the glitch. The timing observations also gave an improved 
dispersion measure of 371.1 ± 1.7 cm^^ pc. 

> 

Q\ ', Subject headings: pulsars: individual (PSR J1718— 3718) 

in 

O I 1. Introduction 



PSR J1 7 18— 3 718 was discovered in the Parkes multibeam pulsar survey (PMPS) by 



Hobbs et al.l ( 120041 ): it is distinguished by its long period, P ~ 3.378 s and very high period 



\-i ' derivative, P ~ 1.598 x 10 Assuming magnetic dipole braking, these parameters imply 
a relatively low characteristic age Tc ~ 34, 000 yr and a very high surface dipole magnetic 
field Bg ~ 7.4 X 10^^ G, one of the highest known for radio pulsars, putting it close to the 
magnetars, i.e., anomalous X-ray pulsars (AXPs), and soft gamma-ray repeaters (SGRs), 
on the P — P diagram. In fact. X-ray emission at a similar level to that from quiescent 
A XPs was detected in Chand ra data from a position consistent with that of the pulsar 



by iKaspi fc McLaughlinI (120051 ). More recent Chand ra observations h ave revealed that this 



X-ray source is pulsed at PSR J1718— 3718's period ( IZhu et al.l 120111 ) . confirming the asso- 
ciation. A spectral analysis showed that the X-ray emission is thermal with a black-body 
temperature higher than normal for rotation-powered pulsars of similar age, strengthening 
the arguments that PSR J1718— 3718 is a quiescent magnetar. 
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Radio timing of this pulsar using the Parkes 64-m radio telescope commenced soon after 
the pulsar's discovery in 1999 and continued with some gaps until 2007 September. In order 
to provide conte mporaneous timing data for the Chandra X-ray observations reported by 



Zhu et al.l ( 120111 ). timing observations recommenced in 2009 January. It was immediately 



obvious that the pulsar had suffered a massive glitch since the earlier observations. In this 
paper we report on the characteristics of this glitch]^ We also give an improved dispersion 
measure (DM) for the pulsar and discuss 3.1 GHz (10 cm) pulse polarization results including 
the pulsar's rotation measure (RM). 



2. Timing observations and analysis 



Following the pulsar discovery, timing observations were made using the center beam 
of the 20-cm Multibeam receiv er and the analogue filterbank (AFB) system used in the 
PMPS jManchester et allboOll ). The AFB system has a bandwidth of 288 MHz centered 
at 1374 MHz and observations were typically of 10 - 20 min duration. These observations 
ceased in 2007 February. In 2005 December, several timing observations were made using 
the Parkes digital filterbank system PDFBl with the 10-cm receiver; this receiver has a 
bandwidth of 1024 MHz centered at 3100 MHz. Several observations of pulsar were also 
made with the 20-cm system and PDFBl between 2007 July and 2007 September. Regular 
10-cm timing observations, mainly using PDFB4, commenced in 2009 January. This system 
gives 1024 bins across the profile and 1024 channels across the 1024 MHz bandwidth. The 
20-cm PDFBl observations were of 10 - 15 min duration and those at 10 cm of 10 min 
duration. All PDFB observations were preceded by a 1-min observation of a pulsed cali- 
bration signa l injected into the feed w hich enables calibration of the system gain and phase 



van Straten et al.ll2010l ). Observations of Hydra A were used to place the data 



response (see 

on a flux density scale. Data were processed using the PSRCHIVE pulsar an alysis system 



2006 



( Hotan et al.l 120041) and th e timing analysis was done using Temp o2 package (jHobbs et al. 



Edwards et al. 2006 ). The Solar-System ephemeris DE405 (Standish 19981 ) was used 



to convert pulse times of arrival (ToAs) at the Observatory to the Solar-System barycenter. 
There were no significant offsets between the ToAs from the different observing systems. 
Results are quoted in the barycentric dynamical time (TDB) system referenced to TT(TAI) 
and quoted uncertainties are one standard deviation. In this paper, we discuss data up to 



^ Given the large gap between the last pre-glitch observation and the first post-glitch observation, we 
cannot rule out that there was more than one glitch. However we consider this very unlikely as no other 
glitches have been o bserved over the 12-y ear data span and large glitches do not normally occur in quick 
succession (see, e.g.. lEspinoza et al.ll201l[ ). 
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2011 March 28. 

Figure [T] shows the variations in pulse frequency v = P^^ and frequency derivative z> 
through the observed data span, obtained from independent fits to short sections of data. 
The top panel clearly shows the large step in u resulting from the glitch (the approximate 
glitch epoch is indicated by the vertical dashed line) and the middle and lower panels show 
that |z>| increased significantly after the glitch. More interestingly, there appears to be a 
significant downward trend in the z> values after the glitch, that is, an increase in 

The changes in u and its derivatives can be quantified by timing fits to the pre-glitch and 
post-glitch data spans; results from such fits are given in Table [TJ Uncertainties (la) in the 
last quoted digit are given in parentheses. This pulsar exhibits a substantial amount of timing 
noise and, for both the pre-glitch and post-glitch data, the timing parameters were obtained 
after whitening of the post-fit residuals by fitting of high-order frequency derivatives as 
indicated in Table [1] A timing position was determined from the whitened pre-glitch data: 
RA(J2000) 17^ IS'" 09!5(6), -37° 18' 49'.'7(23). Thi s position i s con sistent with but less 



precise than the Chandra X-ray position reported by IZhu et al.l (120111 ): we therefore adopt 
the latter for the subsequent timing analyses. The changes in the pulse frequency and its 
derivatives at the time of the glitch (Az/^, Aug and Ai>g) can be estimated from the pre- and 
post-glitch solutions. It is also possible to fit directly for them in a timing solution for the 
whole data span. Rather than deal with the tim ing noise by fittin g of frequency derivatives. 



we have chosen to use the "Cholesky method" ( IColes et al.ll201l[ ) which is implemented in 
Tempo2. This method explicity includes a model for the "red" timing noise and gives more 
realistic estimates of the parameter uncertainties. For PSR J1718— 3718, a fit to the power 
spectrum of the timing residuals showed that it is well-modelled by a power-law with 
a cut-off at 0.2 cycles yr~^ and this was assumed in the Cholesky fit. Results from this fit 
are given in the last column of Table [U A glitch epoch midway between the last pre-glitch 
and the first post-glitch ToA was assumed and its quoted uncertainty covers the gap in the 
data. Post-fit timing residuals are shown in Figure |2l The resulting values of Az/g, Aug and 
Ai)g are consistent with the differences in z/, z> and i) between the post-glitch and pre-glitch 
solutions given in Table [1] when referred to the assumed glitch epoch. Figure |2] shows that 
the next term in the Taylor series, z/', appears to have changed sign from negative before 
the glitch to positive afterward. Further observations are required to verify that this is a 
glitch-related change. 

Comparison of ToAs from 20-cm AFB data and 10 cm PDFBl data recorded between 
2005 June and 2006 February allowed a more accurate determination of the pulsar's DM. To 
compensate for the effects of interstellar scattering, zero phase of the 20cm template profile 
was placed on the leading edge of the profile at about 60% of the peak amplitude, leading the 
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peak by 0.005 in phase. The DM value resulting from the Tempo2 fit is given in Table [2l 



3. Polarization observations and analysis 



All post-glitch 10-cm PDFB4 observations recorded full Stokes parameter data. A total 
of 19 such observations were summed to give a total integration time of approximately 3.1 h. 
After calibration for instrumental gain and phase, the observations were summed in time and 
then summed in frequency with a range of rotation measures (RMs) from —2000 rad m~^ 
to +2000 rad m~^ in steps of 40 rad m~^. A clear peak in the linearly polarized intensity 
L = (Q^ + U'^Y^'^ was observed at about —155 rad m~^. Using this RM, the data were 
summed to form profiles for two frequency bands. The weighted mean position angle (PA) 
difference between the two bands across the pulse profile and a corresponding ARM were 
then computed. This process was iterated until convergence to give the final RM value given 
in Table H 

The final 10-cm polarization profiles shown in Figure [3] were then formed by summing 
in frequency using the final RM and DM. Table [2] gives the mean flux density of the pulsed 
emission, the pulse width at 50% of the peak intensity and the fractional polarizations, where 
the represents a n iean value over the pulse profile. Noise biases in (L) and {\V\) have 



been subtracted (see lYan et al.ll201ll ). T here is no evidence for any change in the 10-cm 



pulse profile at the time of the glitch (cf., IWeltevrede et al.ll201ll ) but, because of the low 
S/N ratio, a small change cannot be ruled out. 



Discussion 



The fractional pulse frequency change for this glitch, Aug/u = (33.25 ± 0.01) x 10 ^, is 
the largest ever observed for any pulsar, magnetars included. The largest previo usly known 
was fo r PSR B2334+61 which had a glitch with Aug/u ~ 20.5 x 10"^ in 2005 JVuan et al. 



2010al ) and for magnetars, the largest glit ch so far observe d is for IE 1048.1—5937 (PSR 
J1048-5937) with Aug/u ~ 16.3 x 10"^ jDib et all 120091 ) 1^ As Figure □ and the glitch 
parameters in Table [T] show, there was a step change in both z> and i) at the time of the 
glitch. The observed value of Aug/ v ~ (6.0 ± 0.7) x 10"'^ is a little less than values typically 



ilsrael et al.l (|2007l ) claimed a very large glitch (At/„/i/ ~ 65x10-^) i n the AXP CXOU J164710. 2-455216 



(PSR J1647— 4552) but this claim has been disputed bv lWoods et all (j201lh who saw no evidence for such 
a glitch in a more complete data set. 
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observed in other large glit ches, but certainly within the observed range (e.g., lYuan et al. 
2010bl : lEspinoza et aLlboilj ). 



Typically, after a large glitch, the spin-down rate relaxes back toward the pre-glitch 
value z>o approximately following the relation 



-t/Td 



where Aup is a quasi-permanent change to z> at the time of the glitch a nd Au^ is the am- 
plitude of a compon ent which decays away on a timescale (see, e.g., lYuan et al.ll2010al : 
Espinoza et al.l l201l[ ) . In some cases the Ai^p is not really permanent but slo wly decays in 
an approximately linear fashion (e.g., IShemar fc Lynelll996l : I Yuan et al.l 1201 Obi ) and in other 



cases there is little or no exponential decay. If present, the exponential term implies an 
increase in i) at the time of the glitch which decays with the same time constant. 

The post-glitch behaviour of PSR J1718— 3718 was very different to the types of relax- 
ation described above, all of which imply both Ai)g > and a post-glitch value of i) > 0. 
For PSR J1718— 3718 there was a significant decrease in i) at the time of the glitch and no 
evidence for any exponential recovery. The step decrease in i) at the time of the glitch is 
shown by the post-glitch downward trend in z> seen in Figure [H by the difference in the fitted 
i) values given in Table [1] for the pre- and post-glitch data spans and by the value of Ai)g 
from the glitch fit. It is possible that there was an exponential decay in z> that we missed 
because of the large data gap around the time of the glitch but, since we see no evidence 
for it in the post-glitch data, its decay time constant would have to be < 250 d. The change 
in i) has persisted for more than 700 days and seems quite distinct from the period fiuctua- 
tions observed prior t o the glitch (Figure [1]). None of the 32 large glitch events discussed by 
Espinoza et al.l (120111 ) show this type of post-glitch behaviour. 



Glitches in AXPs tend to have unusual proper ties compared to those in normal radio 



pulsars (e.g.. lDib et al.ll2008t iLivingstone et al.ll201l[ ). With its relatively long period, strong 



implied surface dipole magnetic fi eld and X-ray emission, PSR J1718— 3718 ha s been identi- 
fied as a possible quiescent AXP (jKaspi fc McLaughlinI l2005l : IZhu et al.ll201ll ). The unique 
post-glitch timing behaviour of this pulsar may be yet another example of these unusual 
properties. 

The recoveries observed after most large glitches are attributed to the changing torque 
on the neutron-star crust as differential rotat ions of the internal superfiuids and the crus t 
relax back toward their equilibrium values (e.g. lAlpar et al.lll993l : lAndersson fc Comerll2006l ) . 
However in the case of PSR J1718— 3718 we do not have a relaxation, we have the opposite 
— the effective braking torque increased by a small amount at the time of the glitch and 
has continued to increase since then. The continued increase implies a negative braking 
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index — the post-glitch effective braking index is about —150 whereas before the ghtch it 
was approximately —15. The characteristic timescale for the increase in torque is relatively 
short: z>/i> ~ 450 yr. 

A secular increase in the post-glitch braking torque is not predicted by standard glitch 
models. In the superfluid model it implies that the differential rotation of the crust and 
superfluid interior, decreased by the glitch, continues to decrease despite the secular slowing 
of the crustal rotation rate. The short characteristic timescale for the post-glitch increase in 
braking (or equivalently, the large negative value of the post-glitch braking index) implies 
that the effect is not purely electromag r ietic . If the glitch resulted in, or was caused by, a 



crustal shift as envisaged by 
in the interior of the star (IRuderman et al 



RudermanI (1991 ) , then it is possible that complex interactions 



19981 ) could result in an increasing braking 



torque. Another possibility is that a heat pulse associated with the ghtch could result in 
an increasing particle flow in the magnetos phere as the surface t emperatur e increases . Such 
thermal flows can have a timescale of years (ILarson fc Linkll2002l ). However IZhu et al.l (120111 ) 
found no evidence for any incre ase in the X-ray flux associated with the glit c h. Recent s tudies 



of neutron-star superfluidity (jPeralta et al.l l2006l : lAndersson et al.l 120071 : iLinkI 120111 ) have 



shown that the onset of turbulence in the interior superfluid can modify the spin behaviour 
of pulsars following a glitch. Specifically, the tangled vortices result in an increase in the 
frictional coupling between the superfluid and the rest of the star and hence a decrease in 
the relative velocity of the two components. If such turbulence were induced by this large 
glitch, then it is possible that the enhanced dissipation could further decrease the differential 
velocity resulting in an increasing effective braking torque. 



We thank our referee for pointing out the possible role of superfluid turbulence in mod- 
ifying the post-glitch recovery. RNM is a CSIRO Fellow based at CSIRO Astronomy and 
Space Science. GH is supported by an Australian Research Council QEII Fellowship (project 
#DP0878388). We thank our coUeages for help with the Parkes observations of this pulsar. 
The Parkes radio telescope is part of the Australia Telescope, which is funded by the Com- 
monwealth of Australia for operation as a National Facility managed by the Commonwealth 
Scientific and Industrial Research Organisation. 
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Table 1. Timing parameters for PSR J1718— 3718 





Pre-glitch 


Post-glitch 


All data 


"DA { To^^^^^a 
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Nr of frequency derivatives 


8 


3 


2 


Epoch of pulse freq. (MJD) 
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55250 
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Data span (MJD) 
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Rms timing residual (ms) 
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Reduced x^/d.o.f. 


2.3/83 


4.3/26 




Glitch epoch (MJD) 






54610(244) 


Az/, (Hz) 






9.842(3) X 10-6 


Avg (10-15 Hz2) 






-0.84(9) 


Ai>3 (10-24 Hz3) 
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'^Position from Chandra X-ray data ( IZhu et al.ll201ll ) 



Table 2. Interstellar and profile parameters for PSR J1718— 3718 



Parameter 



Value 



Dispersion measure (cm ^ pc) 

Rotation Measure (rad m^) 

Mean flux density at 3100 MHz (/) (mJy) 

Pulse width at 50% of peak (ms) 

Fractional linear polarization (L) / (J) (%) 

Fractional circular polarization {V) / {!) (%) 

Fractional abs. circular polarization (|V^|)/(/) (%) 



371.1(17) 
-160(22) 
0.30 
41 
29 
5 
7 
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Fig. 1. — Variations of pulse frequency v and frequency derivative i> for PSR J1718— 3718. 
In the middle plot, the mean post-glitch frequency has been subtracted to show more detail. 
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52Q00 540Q0 56000 
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Fig. 2. — Timing residuals for PSR J1718— 3718 after fitting for pulse frequency v and its 
first and second time derivatives and for the jumps in these quantities at the time of the 
glitch using the Cholesky method. The glitch occurred in the gap between the higher-quality 
PDFB4 data and the earher AFB data at about MJD 54600. 
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Pulse Phase 



Fig. 3.— Polarization profiles at 10 cm (3100 MHz) for PSR J1718-3718. In the lower part 
of the figure the heavier full hne is the total intensity (Stokes /), the hghter full line is the 

linearly polarized intensity (L = {Q^ + f/^)^/^) and the dotted line is the circularly polarized 
intensity (Stokes V). The upper part of the figure shows the position angle of the linearly 
polarized part {if) = 0.5 tan~^(t//(5)). 



